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I - AERODYNAMIC THEORY OF OSCILLATING' WINGS OF FINITE SPAN 


By Eric Reissner 


SUMMARY 


A formula is derived for the pressure distribution on an oscillat- 
ing lifting surface of finite span under the assumption of a ratio of- 
span to average chord (aspect ratio) that is not too small. The range 
of validity of this formula, so far as aspect— ratio limitations are con- 
cerned. is not less than the range of validity of lifting— line theory 
for the non-oscillating wing. 

It Is found that the effect of three-dimensionality of the flow may 
be incorporated in the results of the two-dimensional theory by adding a 
correction factor a to the basic function G’(k) of the tvor-dimensional 
theory. 

The correction term a is a function that depends on wing plan 
form, wing deflection function, and reduced frequency k. Its determi- 
nation requires the solution of an integral equation which is similar to 
the integral equation of lifting-line theory. 

The present report concludes with an explicit statement of the form 
which the results of the theory assume for the spanwise variation:., of 
lift, total moment, and hinge moments on a wing which is oscillating in 
bending, torsion, and aileron and tab deflection. 

Methods for the numerical evaluation of the results obtained as 
well as numerical applications to specific problems are given in part II 
of this re port. G 


INTRODUCTION 


The present report deals with the linear aerodynamic theory of 
oscillating airfoils of finite span. It contains the outcome of attempts 
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to obtain simplifications and extensions of earlier results given in ref- 
erences 1 to 3. The guiding thought in the developments is to deduce, 
from a rigorous formulation of the problem of nearly plane lifting sur- 
faces. simplified results in applicable form which depend on the assump- 
tion of sufficiently large aspect ratio. It is found that such results 
can be obtained and that their range of applicability, so far as aspect- 
ratio limitations are concerned, is at least as inclusive as that of 
lifting-line theory for the airfoil in uniform motion. 

In reference 1 wings of rectangular plan form were considered and 
expressions were obtained for the spanwise distribution of lift and mo— 
ment on the airfoil for arbitrary, small periodic displacements of the 
points of the wing surface. In reference 2 these results were extended 
to wings of arbitrary plan form. In reference 3 expressions were ob- 
tained, on the basis of results in references 1 and 2, for the spanwise 
variation of aileron hinge moment and tab hinge moment. Moreover, it 
was established that the effect of finite span could be incorporated 
into the results of the two-dimensional theory, as given in references 
h and by a modification of the basic function C = F + iG in an 
explicitly specified manner. The term to be added to the function C 
in order tc account for the three-dimensionality of the problem was 
found to be the same for lift, moment, and hinge moments but to depend 
on the nature of the motion being dealt with. Determination of these 
three-dimensional correction terms requires the solution of an integral 
equation for the variation of the circulation along the span. The deri- 
vation of this integral equation is an important part of the work. 

In this report the foregoing results are obtained in what appears 
at present to be the simplest possible way. In- particular, a consider- 
able reduction in the necessary analysis is accomplised by showing that 
the effect of i'inite span manifests itself in the expression for the 
chordwi3e pressure distribution solely by a modification of the function 
C. With this result, use may be made of the known formulas for lift and 
moments of the two-dimensional theory in order to establish the final 
results of the present theory without further integrations. 


SYMBOLS 


U velocity of flight 

x, y, z Cartesian coordinates 


u+U, v, w components of fluid velocity 


t 


time 
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P 

p 

0 

*i(y) 
*t (y-) 


b 

b. 


w. 


a 


V R 


CO 


r 

n 

K 


y*, t \* 


z > t 


pressure 

density 

velocity potential 

discontinuities of u and v in the xy-plane 
coordinate of leading edge 
coordinate of trailing edge 
semichord 

semichord of midspan., 

ratio of span to chord at midspan 

velocity component w at the airfoil 

instantaneous deflection of points of wing surface measured 
normal to xy— plane 

regions in xy— plane (airfoil, wake, remaining) 
circular frequency 

barring of terms, defined by equation (12) 

circulation per unit of span 

circulation function defined by equation (15) 

function defined by equations (39&/b) and occurring as 
kernel of the integral equation (62) 

dimensionless spanwise coordinates; y* = y/b Q , rf* - q/b 0 

dimensionless chordwise coordinates defined by equation 
(40b) 


Zja dimensionless coordinate of midchord line; z m =(xj + x^) / 2b Q 

k reduced frequency; k = cn b/U 

k Q reduced frequency at midspan 
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Q 

Aj ) Ag 

b 

Q (a '> 


s 

C(k) 

0 

a 

c 

e 

d 

f 

1, m 

h 

a 


P, 7 


L 


Mp, ^ 

F 


^0 ^m 

function defined by equation (Vf) 

functions defined by equations ( 52 b) ana ( 53b) 

function defined by equation (59) 

two-dimensional circulation function given by equation (60) 

variable of integration' 

function defined by equation ( 72 ) 

function defined by Theodorsen 

correction term defined by equation (77) 

location of elastic axis in units of semi chord b 

location of aileron leading edge in units of b 

location of aileron hinge line in units of b 

location of tab leading edge in units of b 

location of tab hinge line in units of b 

aileron and tab overhang in units of b; l = e — c, m=f— d 

bending deflection of wings 

angle of attack of wings 

aileron and tab deflection angles 

lift per unit of span 

moment about elastic axis per unit of span 
aileron and tab hinge moments per unit of span 
function defined by equation (86) 


I 
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FORMULATION OF TEE PROBLEM 


The wing considered is of arbitrarily given plan form and is placed 
in the path of a uniform incompressible air stream of velocity U. Small 
camber, thickness, angle of attack, and deformations are assumed for the 
wing which, for this reason, may be represented by a plane surface of 
discontinuity of pressure and tangential velocity in the fluid, parallel 
to the direction of U. In addition to the surface of discontinuity 
representing the wing itself, there is admitted the possibility of a 
surface of discontinuity of tangential velocity (but not of pressure) 
which extends downstream from the trailing edge of the wing and. which is 
also taken to be parallel to U. 


Let x and y be the coordinates in the plane of the surfaces of 
discontinuity, the direction of U determining the x-axis, and let z 
be the coordinate perpendicular to x and y. Let u + U, v, and. w 
be the components of fluid velocity in the x, y, and z directions, 
respectively. Let p designate the pressure in tho fluid and p the 
density. With the higher order terms in the velocities neglected, the 
equations of fluid motion and of continuity become 


(a) 


du y du 

1 


dt dx 

P 

dx 


(b) 


dy y ov 1 

dt dx p 



( 1 ) 


dw y dw 1 dp 

dt dx p dz 


du dy dw _ 

dx + dy dz 


( 2 ) 


, Exterior to the surfaces of discontinuity the flow is without 
vorticity so that in termB of a velocity potential 0, 
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(a) 

«.M 

dx 

0 >) 

V = M 

• by 

(c) 

, M 

* 3z 


From equations (3) and (l) an appropriate form of the Bernoulli equation 
follows for the pressure p , 


Let 



■ _ £ . » M + 

uM 



p ot 

dx 


(a) 

7 = (Ziu) ~ 1 

t-g 

\ 

Jz-o 

(V) 

5 = (Av ) z _ 0 - | 

[ A 7- 

V oy 

-^z=o 


(4) 


(5) 


stand for the discontinuity distributions of the velocity components u 
and v ever the surfaces of discontinuity in the xy— plane . According 
to the Biot-Savart theorem the velocity component w in the interior of 
the fluid is then given by the integral 



y_i | > T)j t)(x - I ) + s( \ , Tb t)(y - q) d ^ 

Ux-t f* (y -->) 2 + z s 'Y /z 

^ J 


( 6 ) 


From equations (5) it follows that 7 and 5 are related to each 
other in the form 


by _ b 6 
by bx 


( 7 ) 
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From equations (4) and (5)' it follows that the discontinuity of 
pressure Ap is given by 



( 8 ) 


The problem of the oscillating lifting surface may now be formu- 
lated. Designate the coordinates of the leading and trailing edge of 
the wing by xj( y) and x.j. (y), respectively, with Xj(0) = — b Q and 

x t (°) = b^ so that b G represents the semichord at midspan. Let sb Q 

stand for the length of the semi span. Distinguish the following three 
regions in the xy— plane; (l) the airfoil region R a , (2) the wake 

region which is the strip of width 2sb 0 extending downstream from 

the trailing edge, (3) the remaining region R r . Then there is the 

condition of continuous flow outside the airfoil and wake regions, in R r , 

7=5=0 (9) 

and the condition of continuous pressure in the wake region R w , 

Ap = 0 (10) 

At the airfoil the condition of tangential flow is tc be satisfied. If 
Z a (x, y, t) stands for the instantaneous deflection of points cf the wing 

surface measured normal to the xy— plane, this condition in the region 
R a is 


w a = w(x, y, o, t) = 


c&a 

5 t 




+ U 


a 
d x 


(id 


By introducing equations (9) to (ll) in the integral representation 
for w, equation (6). and taking account of equations (7) and (8), there 
may be obtained the general integral equation of lifting— surface .theory 
for the distribution of the velocity discontinuities 7 a and S a over 

the airfoil region. With the solution of this integral equation, 
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equation (8) furnishes the chordwise pressure distribution A p a . Appro- 
priate integration of A p a results in expressions for spanwise lift, 
moment, and aileron and tab hinge-moment distributions. 

Further development is here based on the assumption of simply 
harmonic motion. 


z a u, = z a (*> y)e iU)t 


( 12 ) 


As a consequence of the linearity of the theory, it follovs also that 
velocities and pressure become products of e 1£U ^ and of amplitude func- 
tions independent of time which are designated by bars. 

The first step in the deduction of the final form of the integral 
equation of the problem consists in the determination of the values of 
7 and 6 in the wake region, by means of equations (8), (7), (9), and 
(10). Equation (8) takes on the form 

x t X 

0 = its ^ 7 a dx' + j 7 W dx' ^ + U 7 W ( 13 ) 

: ; * 1 . *t 

If the circulation F is defined by 


x t 

F = P 7 a dx* 

x z 


and the circulation function Q by 


( 14 ) 


Q 



. Oi 

1 — x+ 
e U t 


( 15 ) 


it follows from equation (8) and its differentiated form that 7 W is 
given by 
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7 W — 5 - b d 


0) --- -i 
- f 2 e 
U 


<n 

U 


x 


(16) 


When 7v of equation (16) is substituted in equation (7), there is ob- 
tained for 5 W , 


6 - [' f d*' 

w j 5y 


—CO 


= b. 


— .00 
dft “i r; x 

e u 

dy 


(17) 


Equations (17), (l6) , and (ll) are now introduced in the integral 
representation for w, equation (6), in which the coordinate z is 
made to approach zero. As was shown in reference 1, the two limit proc- 
esses of integration and of letting z approach zero may be interchanged 
if the Cauchy principal value of the integrals is taken at the points 
{ = x, q = y for which the integrand becomes infinite. There results 


v a (x, y) = 


-1 TP ?a( £ > ti)(x - |) + & a ( i , rjHy-rO 


•Jj 

*a 


i (x - i) s + (y -q) 


2 j 3 / 2 


d i dp 


_^o 

k n 


.CD& , (D — / \ > IK 4Q , 

f'p -1 TT -i 7 t (q)(x - |) + gjy- (y - tj ) 


JJ 

R. 


w 



(x - i ) 2 + (y - p) 2 

3 / 2 

1 

L 



d | dp (18) 


Equation (l8) holds for all points x, y inside R a . Its left-hand 
side is the given function 


Sz g 

w a = icD Z a + U — 


(19) 


and, according to equations (7) and (9), & is expressed in terms of 

cl 

7a by means of the formula, 
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( 20 ) 


In addition, there is imposed the Kutta-J oukowsky condition stating 
that 7 and S are finite along the trailing edge of the wing. The 
amplitude of the chordyise pressure distribution A p a in terms of the 
function 7 a , for which equation (18) is to be solved, is given by 


A p 
*a 

P 


x 



7 a cbc' + U 7 a 




( 21 ) 


The general problem is the solution of equations (l8) to (20) for 
arbitrary shape (plan form) of the airfoil region R a . Results equiva- 
lent to this solution are known for the following two special cases: (l) 
the two-dimensional problem to which the general problem reduces when 
the region R & is the strip jxj < b = b 0 and w a of equation ( 19 ) 

is independent of y (see for instance reference 5), and (2) the problem 
of the wing with circular plan form (references 6 and 7). 

The purpose of the remainder of this report is to obtain approximate 
solutions of the problem which are applicable subject only to the re- 
striction that the airfoil region R a is of sufficiently large aspect 

ratio, a restriction that is roughly equivalent to requiring that the 
ratio s of span to chord at midspan is sufficiently large. The nature 
of these results is, presumably, that of the dominant terms of an asymp- 
totic development of the exact solution in terms of the parameter s. 
However, no investigation of this aspect of the problem is made in the 
present report. 

Concerning the range of validity of the results as far as aspect- 
ratio limitations are concerned the statement may be made, in view of 
the nature of the analysis, that the results obtained in this report for 
oscillating wings are certainly applicable for all wings for which 
lifting-line theory is considered applicable in the case of uniform mo- 
tion. There is some reason to believe that the results for oscillating 
wings may have a somewhat wider range of validity than lifting— line 
theory, as the chordwiee waves which occur in the vibration problem may 
be responsible for a reduction of the effective chord length while at 
the same time having no influence on the effective span length. 
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DERIVATION OF SIMPLIFIED INTEGRAL EQUATION OF LIFTING-SURFACE THEORY 
Equation (l8) may be written in the form 


x t> - S V ^ y ^ eb Q ; w a = I x + Ij 


( 22 ) 


where 


^0 x t^L 


(a) lx = 


k n 


(b) i 2 = -2 


f 

f 

J , 

-8b 0 

x z(ri) 

m 

o' 

0 

OO 

r 

f i 
/ 

.1 . 
-eb o 

/ 

x t (n) 


(x -I ) 2 + (y - q) ? 


3/2 


aUt) 


(23) 


i - t -i ^ Q (n)(x - i ) + Q (t))(y - q) 


e U 


(x - i ) 2 + (y - r|) 2 


3/2 


a bd'n 


and where ; from equation (20), 


§ a ( i, t) = 


zjygCVj, n) 

d n 


a 


x zU) 


The decisive step in the present solution of the problem as ex- 
pressed by equations (22) and (23) is the reduction of the two- 
dimensional integral equation to 2 one-dimen3ional integral equations 
which can be solved successively. For this reduction it is assumed that 
the airfoil region is sufficiently elongated in span direction, that the 
rate of taper is moderate, and that the velocities do not vary appre- 
ciably along the span over distances which are a fraction of the mean 
chord. With these assumptions it may be postulated that at every span- 
wise section the contribution to the normal velocity w a induced by the 
tangential velocities at the airfoil is approximately as if every 
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spanwise section were part of a. wing without taper and as if the flow 
were two-dimensional. Equation (23a) then becomes 


I, — 


Ll & 


t 7 a ( I , y)(x - i )d t dr) 


4 n 


f f 

-00 xj(y) 


(x - i ) 2 + (y - ■q) 2 


x t (y) _ 

-1 r i » y) 

2 * J 


X 


- i 


*i(y) 


3 /s 


dfc = I. 


(24) 


A corresponding approximation I 2 for the integral I 2 is ob- 
tained by a less, direct procedure . Write 


-b 

Is = “ ( I 


-( 

Jt \ 


. 00 _ 

3 1 - 1 4 

3 U 


) 


(25) 


where 


I = 
3 


sb„, 

O 00 

n n 


./ 


_i ® ? 

x • i> 

e U 


Sllhlll-z j 1_ a | dtl (ae) 


| (X- if * (y-n) 2 1 3/a 


-sb 0 *t(l> > '" *' ' ” J 

sh 


0 00 


. 00 

—1 


J 4 = 


r 


,/ ./ c , 

- sb o x t^ 


H Q inl.il ~ U 


(x - i ) 2 + (y -q) 2 13/3 


d | dr) (27) 


Write next 


sb 0 o° 


h » 


./ j ~f . 

-sb 0 x — sb 0 x 


sb 0 x t (r)) 


(28) 
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In the integral extended from x to x t (p), that is, over part of the 

airfoil region, the same approximation can be made as in the expression 
for Ii. If this is done, its contribution to I vanishes and 

3 



~i 

e U O' (n)(y - p) 


( 


X -t) Z + (7 - r\Y 


3 /a 


-d| dp 


With a new variable of integration X = % - x 


I 


3 



MI 

dp 


r 


O) 

00 ~i £ X 
e 


(y - 1 


x 2 + (y - p) 


2-3/2 


In order to approximate the remaining integral I 4 this integral is 
separated into two parts, one part which equals the value of I 4 in the 
two-dimensional theory and the other part being the difference between 
the two-dimensional and the three-dimensional value of I 4 , 

I 4 = I 4 (s) + ai 4 (30) 


00 03 



CD fe 

e ^ (x - t )d£ dp 
(x - t ) 2 + (y - p ) 2 1 


= 2 n (y) 





(x - t ) 


dt 


(3D 


The factor Q (y) may be written in the alternate form 



14 


NACA TN No. 1194 


ft (y) 





2 



ft (t|) — — drj 

y - tv 


(32) 


if account is taken of the fact that ft ( ± sb 0 ) = 0. 

If equations (31) and (32) are combined, the value of I 4 ^ s ^ may 
be written 


I 


4 


( 2 ) 



-i ® i 


\j > 9. tul ly ~ n 1 

(x - i) (y — TJ ) 


dl drj 


(33) 


The next 3 tep is transf ormation of I 4 by integration by parts 

with respect to t] . Before this step is carried out it is noted that it 
is consistent with the preceding approximations to replace in equation 


: r 

I 

x t (y) 


. CD 

— 1 r; 
e 


<* - £) 


x t (rj) by 

x^(y). Then 

sb n 


[ r _ 

ft (rj) drj 1 

\ / 

r ■ 


v : 

(x - n 2 + (y - ti) 2 

1 3/3 J 


AS 


00 

n 


CD 


(x-n 


x+(y) 


eb. 


ft' (t)) ( y -- T\) dr) 


-sb. 


(x- nVU - U 2 + (y - rj ) 2 -J 




( 34 ) 


as the integrated portion of the inner integral vanishes because 
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Equation (35) for A I 4 may again be separated into two integrals, in 
the sane manner as I 3 of equation (28). If this separation is made, 
the integral between the limits x and x fc (y ) may again be neglected 
with the result that 



and with \ = i — x 
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Now, "by combining equations (36) and ( 31 ), there follows: 




00 
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For the applications it is convenient to transform K(x) by an inte- 
gration by parts of the first integral as follows: 


» -iX, 

1 e 6X ~i\ 

= e 


(X 2 + x 2 ) 


3 / 2 


X? -f~ X 2 


+ c(x) 


+ x c 


+ i / e 


-LX 


( 


x 2 ,J X s + x 2 


) 


+ c ) dX. 


The right tide assumes a definite value only if 


c(x) = — x" 


and then 


00 

r 


-ix, ^ 

e dX 


2 L .-2\3/3 


(X S + X ) 


1 r e- ix 

2 + * / 2 
X / X 

0 


( 


f X 2 + x 2 


-0 


dX, 


By introducing this formula in equation (39a) there follows: 


e~ 1X ( 1 - ZE EZZ z l *! ) IX ( 39 b) 


Combination of equations ( 38 ) , (24), and (22) results in the approxi- 
mate integral equation of lifting-surface theory, which is the basic re- 
sult of the present work. Because the establishment of this equation 
depends on the condition of a ratio of wing span to average wing chord 
that is not too small (so that the wing plan form has the appearance of a 
strip), it might be considered to call this equation the integral equa- 
tion of lifting-strip theory to distinguish it from the general equation 
of lifting-surface theory on the one hand, and from lifting— line theory 


K(x) = 


cub 


_o 1 
U x 


cub 


o 1 
U x 


00 

n 
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for the stationary sufficiently elongated wing oh the other hand. That 
it is not possible to speak of a lifting— line theory for the oscillating 
wing follows from the fact that according to equation ( 38 ) and (22) the 
effect of spanwise variation of circulation is. hot uniform across the 
chord but varies in a manner depending on the frequency co. 

In order to obtain a normalized form of the equations of the prob- 
lem, the following dimensionless variables may be introduced: 


(a) y* = yA 0 

x ~ \ ( x t + »l) 

(b) * - 

| ( x t - Xj) 


(40) 


As the space coordinate z does no longer occur, from equation (18) on, 
introduction of a dimensionless variable z at this stage will give no 
rise to doubts as to its meaning in the subsequent developments. Write 
as abbreviations 

(a) 


0 >) 

The quantity z m is a measure of the local sweep of the midchord line . 

The quantity b stands for the magnitude of the local chord so that 
b(o) = b Q . With these symbols equation (40b) becomes 


~ Ut ~ x l) = 13 


(41) 


2b“ (x * + - *» 


x = b c z m + bz 


(42) 


Also let 


k = 

U 


k -2a 

° ~ u 


^m " ^o z m 


(43) 


Wherd k is the local reduced frequency and k the reduced frequency at 
midsphn. 
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Upon introducing equations (40) to (4-3) in equations (38), (2.4), aid 
(22) there follows as the normalized form of the integral equation to he 
solved 


-w a (z> y*) 


; X V 1 ” y a ( £ >-r*> 
‘ 2 *_/ * - t 


it 


ikjn 

2~w 


n(y*) 





— s 


Equation (44) has first been established in reference 2. For the 
wing of rectangular plan form (k = k Q , k^ = 0) it reduces to an equa- 
tion given in reference 1. - 

The problem from here on is to solve equation (44) fox the function 
7 a , to obtain an equation for the span vise variation of fi , and to 

express the chordvise pressure distribution in terms of w a . 

In terms of the variable z as defined by equation (42) the func— 
tion q of equations (l4) and (15) is to be written as 


q (y*) = ~ e ' L ^ k+kia ^ J 7 a ( it y*)<U ( 45 ) 

-1 

and equation (21) for Ap a becomes 

z 

- = ik / 7 a <iz' + 7 a ( z, y*) (46) 

p U / 

To be established is the manner in which the presence of the term 
containing dQ/dq* in equation (44) modifies the result of the two- 
dimensional theory for Ap a . 
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CALCULATION OF PRESSURE DISTRIBUTION 


The next Btep in the analysis is the determination of y a from 
equation (44) in which the following abbreviation is inserted: 


s 



-s 


dO 

drr 


K ,[k 0 (y* - H*)] 


dn x 


(47) 


Use is made of the following pair of formulas of integral-equation 
theory: 


l 

g(z)=— f d t, f(l) finite 

2 rt / z - i 

— 1 


f(z) 


£ IkzJk / /1J_L «LLL d i 

n y 1 + z J J 1- ( z - ( 


(48) 


It may be noted that a derivation of the results of the two-dimensional 
-theory in this manner has been given in reference 8. 


By applying equations (48) to equation (44) there follows, 
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y =-2 /m \ r 

7a K / 1 + Z L I 


11 ,/TTT w a(^ y*) a , 

/ i - t z~t 


/ 


/i + (, 


j_v — ik * r 

+ fie <t / r 

2 n J / 1 - '•> 

-l • 


— ikX , i- 

0 dX ) 4 ' 


x / z - i 


-1 


— ik„ 

/ 1 - £ z - $> 


i; 11 - f /^±4 zri« 

2 n 


(49) 


The double integral in equation ( 49) is reduced to a single integral as 
follows : 


f /— ( / 

L / i - c v i 


-ikX , „ 

e ^ d 5 


dX 


UA j 

t- X ^ 


z - b 


= /' e" ikX j / /ll 

J L J / l - 


1+5 d 5 


L / 1-5 (5 - x)(z - 5) J 


! dX 


P — ikX j" 1 

f Aj L 

i ( i i "\ 

1 - i 1 /| 

5 i 

1 

-I V — ■- 

c 

i — 
N 
! 

J / 1 - 

T 

5 \ 5 - x z - 5 / 

J 


dX 


-1 


-ikX, 


1 . 


j f /X + 1 \ 

! « V 1 - / 7 J- K 


I Jl \ M / 

- x L x / x 


dX 


(50) 


By introducing equation (50) in equation (49), 
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(5D 


According to equations (46) and (45) it is necessary. to integrate 
equation (51) between the limits -1 and z, with z' as the variable 
of integration. Hy interchanging the order of integration with respect 
to z', £, and X, respectively, the following integrals occur (see 
reference 8) : 





(52a) 


z 



dz 1 


1 + z' X — z* 


— l 


— + sin 1 z + ] — — A (z, X) 

2 J X + 1 8 


(53a) 


The functions are defined by 




1 1-zt ±jJ 1 - t 2 >/ 1 - z 2 

2 n 1 - z t - y'l _ £ 2 fx _ z 2 


(52b) 


A (z,X) 
2 


2 tan” 



(l — z) (X + l) _ 
(1 + z)(X - 1) J 


(53b) 
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The following formulae will also he needed: 



With equations (52a) and (53a) there follows from equation (51), 



-1 


Equation ( 54) leads immediately to the integral equation for the 
spanwise variation of f. . Let z = 1 and there follows, in view of 
equation (^5) and in view of the fact that A 1 (l, t ) = 0 and 
A 2 (1, X.) = -it, 
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Equation (55) can be simplified by means of the following 
of the theory of Bessel functions: 


CO 



1 


~ik 

= - “ ( hJ 2 ) (k) + iH 0 < 2 ) (k) ) - ®— 
2 \ / xk 


1 

(b) j J~~r ® ik ^ d £ - # ( j 0 ( k ) ~ iJ i ( k ) ) 

Hence 


Q = 


0 b ^ k+k m) 
-2 — e 

*0 


1 



- 1 


J , - ik f n 
+ ikQ e l — 
l 2 


Hi 


(2) 


+ i H ( 


(2) 


-ik 


ik 


1 b - ik 
> Qe 


J b 


(Jo - 


(55) 

relations 


(56) 


Ui) (57) 
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and by canceling two of the Q-terms against eaclj other 

ik m / / i-li _ f 

ii y 1- 

hJ 2 ) + in 0 ( 2 ) 


Q - 


J 0 - iJ l 


nik 


H x ( 2 ) + iH 0 ( £ ) 


k' 


2b _ 
h Q r nik 


(58) 


Define now a function jj, by the relation 


H = - 


J Q - iJi 


J 0 “ iJ i 


nik 


+ iH Q ( 2 ' 


itk 


(J 0 - yi) - t (Ji + y 0 ) 


( 59 ) 


and write as an abbreviation 


/i+ i 

n( 2 ) « 4 1_ e ikm -r* y 1 ~ £ 

b 


Wo d t 


Q 


nik | + iHo^ 2 ^ 


(60) 


Equation ( 58 ) then assumes the form 


n (y*) + n(k) —■ 2Q(y*) = 5^ (y») 
*o ■ 


( 61 ) 


Equation (6l) is the integral equation for the spanwiee variation of 

— (2) — 

circulation. The function Q is the distribution of u according 

to the two-dimensional theory. The Q-terrn represents the influence of 
finite span. With equation (47) for Q there may be written 
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s 

£5 (y*) + n(k) f “7 E [ 

b 0 J *n" 


1r ( xr^ 


(62) 


The next step in the analysis is to subtract equation ( 55 ) mul- 
tiplied by -( — + ein 1 zl from equation (5^)' 

Jt \ 2 / 

z 

7 az*— - ( - + sin -1 z) © _i ( k+k m) q 

a rt V 2 / b Q 

-a 




(63) 


The second and third integral in the right-hand side of equation (63) 
are now integrated by parte. 
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7o dz 1 - - ( -2 + sin 1 z — 

Jt V 2 / b 




- 1 


o r ik — -ik- 

- / Aj v a d ^ 2 Qe 


Jt . 


~ 1 


xf 


- + sin 1 z + A (z, X) 
2 3 


— ikk % ® 

1 1 

-ik J 


CO 


ik_ ikm ['• 5A_ ^- ikX 


ft e 


2 e 


5k — ik 


dX 


-ikt 


- . , / —IE 

-W”(a 

TC 1 —IK 


P - 1 e — at ) 


J a £ -ik 

-l - 1 


(64) 


According to equations (52b) and (53b), 


A x (z, ± 1) = 0 A 2 (z, 1) = 0 


a / v / Jt . -1 \ ‘ 

A (z, oo ) = — i — + am z 

a \? / 




(65) 


The third of these relations is verified by means of the identity 


2 tan 




— z it . — i 

sm ’ z 


+ z 2 


By means of equation (65), and by noting that b Q /b = k Q /k, it is 
seen that the integrated terms in equation (64) cancel and 
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7a dZ ' = \ J A i *a 

"I -1 


k it 6 


5V a- lkx IX 

dX 


-±- ^a' ik “ d£ 

ik n 2 -7 1 8 £ 


(66) 


Now by combining equations (66) and (51), as in equation (46), there 
follows for the pressure distributions on the airfoil 


Ap 


a 


PU 


2 

it 



i - z A + £ 


i + z y i - t z — ^ 


£-T - ik ^i ) w a d £ 


- 1 




It 


(/Hi/rr 


X + 1 1 8 A a ^ -ikX .. 

s dX 


1 z X ^X 


) 


1 rr _ik i 


+ -p Q e 


•m 


(/ 


T~ /I J_L _1 ^7 \ ~ ik ^ 

1 + a y 1 - t z - £ d£ / 


d£ 


-1 


( 67 ) 


Equation (67) simplifies considerably when equations (52c) and 
(53c) are taken into account. There follows 
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l 



Equation (68) contains the noteworthy result that the term with Q and 
also the term with Q lead to chcrdwise pre ssure distributions which 
are both proportional to J (l - z)/(l + z). Note that according to 
equation (44) these two pressure terms are caused by downward velocities 

00 

across the wing of the form / e —i ^ ( z _ £ ) 1 d £ and e — ' . The 

’“I 

fact that these two different velocity distributions lead to the same 
simple pressure distribution is here arrived at by an analysis which 
does not reveal the inner reason for this occurrence. A modification of 
the analysis so as to clarify this point appears to be worth while, par- 
ticularly as the fact itself accounts for the relatively simple form in 
which the aerodynamic span effect modifies the expressions for the air 
forces and moments of the two-dimensional theory. 

In order to obtain the final form of the expression for A pa, use 
the following known formulas of the theory of Bessel functions: 


(a) 

(b) 


00 



(69) 
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and express the quantity Q in terms of Q by means of equation (6l) , 


—(2) - 


2u (lc) k/k 0 


(70) 


Then write 


£*a 

PU 


■ ! / (/it! /rtf ttt - “*> ) % a£ * yffi 3 


where 


(71) 


8 ■ ( r B o ( * > fi ♦ ; 5 ') . 


= — e ik “ I" + 

2 0 








(72) 


—(2) 

Introduce the value of Q from equation (60) and the value of |i 
frem equation (59). Then 
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i 



Bjj*) + 1 H 0 ( £ ) 



Introduce in conformity with reference 4 the function 


C(k) = 


H 




Hi^'ck) + iH 0 Cs) (k) 


(74) 


Equation (73) then becomes 

i^-w a dt!C-l + 

“ 1 



( 


n 


Q 


r^j 


- 1 c + 


i Jx y 
j 0 - iJil 


(75) 
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By introducing equation (75) in equation (71) there follows: 



Equation ( 7 6) is the general result to be established. A corre- 
sponding formula for the two-dimensional theory is given in reference 8. 
Equation (76) shows that the aerodynamic span effect in the present 
theory, manifests itself solely by a modification of the basic function 
C(k) . An additive correction term occurs, which is given in the form 


a 




r c ( k ) + 

L 


iJi(k) ' 
J 0 (k) - iJi(k) _ 


(77) 


and which is seen to depend on the ratio of three-dimensional to two- 
dimensional circulation function f2. In order to evaluate the correc- 
tion term for a given wing deflection function Z a , it is then 
necessary only to solve the integral equation (62) for the spanwise 
variation of circulation. 
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LIFT AND MOMENTS FOR BENDING, TORSIONAL, AILERON, AND TAB DEFLECTION 

For the applications (see reference 10) the following deflection 
functions are of particular interest: 

Bending deflection: Z a = h(y) = hyf^Cy*) 

Torsional deflection: Z a = a(y)(x - ah) = a^,(x - ab)f a (y*) 

Aileron deflection: Z a = f3(y)(x — eb) for cb < x < b 

Tab deflection: . Z = 7(y)(x - fb) for db < x < b 

£L ^ 

Lift and moment functions according to the two-dimensional theory may be 
found in reference 5. Modified so as to include the '-aerodynamic span 
effect these functions can be written in the following form, as was 
shown in reference 3 : 
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Ma 


2pU 2 b 2 






+ 



2pU 2 b 2 


+ 


+ 


+ 


+ 


= it f — - ( - + a ) ik [c + a \ - 

l 2 \ 2 J |_ h j jb 


rtf — A - — 

1 2 02 2 


V -(;♦*) [ 1 + (;-*) “] [ c + »« 


' I . 

i 


f S 

■ 2 


>3 + m s 2 - ■‘‘V] - (5 + a ) f El + f H 


0 + °p ] } p 


r * 

1 2 


A + ikA - k 2 A_, 
7 3 ^2 y 1 


(§*0 


e x ( 4 ) + ~ Ma ) 


C + a. 


] 7 


( 79 ) 


f- — % + - E ik 

! o “1 O 3 


C + a 


h 


h 

b 


c r 

5 ikB - k 2 B 
! 2 as a* 


+ - E 

1+ f - - a ^ ik 


C + 0 

2 3 

V 2 ' / 

| 

a 


) 5 


fi 

1.2 


B P 3 + ^Bp 2 - k 2 B^ 


+ — E, 
2« “ 


E x + — E 2 
2 s 


C + ap 




B + ikB v - k 2 B 
y 3 7s /x 


2 n 


E, 


E x (d) + E 2 (d) 


C + a 


7 J 


} 7 


(80) 
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2pU 2 b c 


* +il.,(4)lkfc + >J 1 


2 "l 2 


) J h 


+ fa [ikDf. - > 2 B_ 1 + - E (cl) ' 1 + - a Y ikl f C + o_l ) cfc 

\ 2 L “e a i 2 3 \ 2 ^_l a I J 


(— f Db - + ikDp, - k 2 !)^-] +— E (d) E x + • E 2 1 T C- + op 1 ) p 

! ,2 L ' . P2 'YJ 2* ? L 2 _ L _ j 


+ ! a D + ikD ~k% * + -- E 3 (d) E x (d) + — E 2 (d) 

‘.2 L 7 .3 7 2 'i J 2jt L 2 


X C + o^ j 7 


The terms A, B, and D are defined' in reference 5- The terms E are 
of the form 


El - T 10 — ZT 2X Ei(d) - Txo(d) — mT 21 (d) 

E 2 = T X1 — 2ZTio E 2 (d) = Tn(d) — 2mT 10 (d) 

E 3 = T 18 - 2ZT 20 E 3 ( d) = T 12 (d) ~ 2mT 20 (d) 


with the terms T also defined in reference 5. 


terms Oj( j = h, a, 0, y ) as defined by equation (77) are given 
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’j= f c + 


J 00 - iJi(k) 


( ITu) " 1 ) 


nj 


(83) 


The functions Qj^ 2 ^ of the two-dimensional theory as defined by equa- 
tion (60) are given by 


ilc_ _ 
- b 4iCe „ h 

= ~ j — r ik - 

11 t > 0 a 1 ( *>(k) b 




Q ( 2 ) = — ^ iCe 


lk m 


a 




a- 


1 + ( — — a j ik 


a 


r (84) 


ik~. 

q (2) = M -MS®. — I (e x + e 2 - 1 — 

P b o fflx^^k) * V 2 ; 


ikv 

O ( 2 ) b 4iCe 1 Itt< /j\ -m / 1 \ ik 

U„ v ' = — — - 1 E^d) + E 2 (d) — 

2 I 


b o kHi^(k) * L 




The functions Q j of the three-dimensional theory are, according 
to equations (62) and (39b), the solutions of the integral equation, 
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n,(y*) + - wOO 


[J 


s- 

r aSj 


u s y* - 


n t 

-ik / F(k|y* - rj-l )d5 , 

°J |y*-ti*| ol 1 J. 




(85) 


The function u is defined by equation (59) us 


n(kj 


J 0 (k) - i Ji(k) 

nk | [ J Q (k) -yx(k)3 - i [Ji(k) + y 0 (k)]J 


The function F is, according to equation ( 39b ) , 


F(x) 



y~x g + x 2 - 

X X 


) 


dX 


( 86 ) 


This function occurred previously in reference 9> where a different 
theory of the problem of the oscillating wing of finite span was put 
f orward. A discussion of the theory of reference 9 is given in ref- 
erence 1. Tables 1 and 2 contain values of the two functions n and F 
for a significant range of the variables k and x. 

It is apparent that the main task in obtaining three-dimensional 
corrections to the results of the two-dimensional theory consists in 
solving the integral equation for Q . The second part of this report 
(reference 10), which deals with applications of the theory, contains a 
practical method for doing this. 


Massachusetts Institute of Technology, 
Cambridge, Mass., December 3, 1945, 
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TABLE I.- VALUES OF THE FUNCTION F OF EQUATION (86) 


X 

F(x) 

X 

i 

r(y. 

:) 

0.00 

00 — 

1.5711 

2.2 

0.113 

— 

0.425i 

.05 

2.750 - 

1.4681 

2.4 

.097 

— 

.3951 

.10 

2.109 “ 

1.3751 

2.6 

.083 

— 

.3691 

.20 

1.490 - 

1.2481 

2.8 

.072 

— 

.3451 

.30 

1.155 ~ 

1.1461 

3.0 

.063 

— 

.3241 

.40 

.935 - 

1.0601 

3.2 

.055 

— 

.3051 

• 50 

.778 - 

.9871 

3.4 

.048 

— 

.2891 

.60 

.658 - 

.9221 

3.6 

.043 

— 

.2741 

.70 

.564 - 

.8651 

3.8 

.039 

— 

.2601 

.80 

.489 - 

.8131 

4.0 

.035 

— 

.248i 

.90 

.427 - 

.7681 

4.2 

.031 

— 

.2371 

1.00 

.376 - 

.7261 

4.4 

.028 

— 

.226i 

1.10 

.333 - 

.6891 

4.6 

.026 

— 

.2171 

1.20 

.297 - 

.654i 

4.8 

.024 

— 

,2081 

1.30 

.265 - 

.6241 

5.0 

.022 

— 

.2001 

1.40 

.238 - 

.5931 

5.2 

.020 

— 

.1921 

1.50 

.214 - 

.5671 

5.4 

.018 

— 

.1851 

1.60 

.194 - 

.5401 

5.6 

.017 

— 

.1791 

1.70 

.176 - 

.5171 

5.8 

.016 

— 

.1721 

1.80 

.160 - 

. 496i 

6.0 

.015 

— 

.1671 

1.90 

.146 - 

.4751 

CO 

l/( 2 k 2 

) 

i/x 

2.00 

.134 - 

.4581 



TABLE II.- VALUES OF THE FUNCTION 

n OF EQUATION (59) 



k 


|i(k) 

k 

H(k) 


0.000 

0.5000 

— 

O.OOOOi 

0.650 

0.2139 - 

0.06651 

.020 

.4810 

— 

.04231 

.700 

.2062 - 

.06101 

.040 

.4607 

— 

.06651 

.750 

.1991 - 

.05571 

.060 

.4410 

— 

.08291 

.800 

.1924 - 

.05071 

.080 

.4226 


.09421 

.900 

.1801 - 

.04l3i 

.100 

.4051 

— 

.10l8i 

1.000 

.1688 - 

.03291 

.125 

.3857 

— 

.10821 

1.250 

.1436 - 

.01591 

.150 

.3690 

— 

. 1120i 

1.500 

.1218 - 

.00421 

.175 

.3522 

**- 

.1130i 

1.750 

.1027 + 

.00311 

.200 

.3393 

— 

.11391 

2.000 

.0864 + 

,006bi 

.225 

.3268 

— 

.11321 

2.100 

.0807 + 

.0070i 

.250 

.3154 

— 

.11161 

2.150 

.0780 + 

.00721 

.275 

.3049 

— 

.10991 

2.200 

10754 + 

.00721 

.300 

.2955 


.1076i 

2.250 

.0730 + 

.0072i 

.350 

.2787 

— 

.1023i 

2.300 

.0706 + 

.0070i 

.400 

.2644 

— 

.09641 

2.350 

.0684 + 

.0069i 

.450 

.2519 

— 

.09031 

2.400 

. 0663 + 

. 0066i 

.500 

.2408 

— 

.08421 

2.480 

.0632 + 

. 006li 

.550 

.2288 

— 

.07751 

2.540 

.0610 + 

.00571 

.600 

.2220 

— 

.07221 

2.600 

.0591 + 

.00521 



